Quantum correlations in the two-photon decay of few— electron ions 
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A theoretical study of the polarization entanglement of two photons emitted in the decay of 
metastable ionic states is performed within the framework of density matrix theory and second-order 
perturbative approach. Particular attention is paid to relativistic and non-dipole effects that become 
important for medium- and high-Z ions. To analyze these effects, the degree of entanglement is 
evaluated both in the dipole approximation and within the rigorous relativistic theory. Detailed 
calculations are performed for the two-photon 2s 1( / 2 — > lsi/2 transition in hydrogen-like, as well as 



for the ls 1/2 2s 1/2 1 So -> ls 2 1/2 1 So, ls 1/2 2s 1/2 3 Si -> ls 2 /2 1 5 and ls 1/2 2p 1/2 3 P —> ls\ /2 1 S 



ls 2 v 

transitions in helium-like ions. 
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I. INTRODUCTION 

Over the decades, two-photon bound-bound transi- 
tions in atoms and ions have provided a unique testing 
ground for advanced atomic theories. Starting from the 
early work by Gdppert-Mayer [l[ and Breit and Teller 
, a large number of theoretical studies were carried out 
to estimate the total as well as the ene rgy - and anglc- 
diffcrcntial (two-photon) decay rat es |3 - 4lO| . When com- 
pared with experimental data [Ill4l5l |. these studies re- 
vealed important information on relativistic, quantum 
electrodynamics (QED) and many-body phenomena in 
atomic systems. Besides structure-related investigations, 
more recent interest focuses on the quantum correlations 
between the emitted photons, which can be used to probe 
fundamental aspects of modern quantum theory. In a se- 
ries of studies, for example, photon-photon polarization 
correlations were employed to test the Bell inequality 
I n particular, these investigations demonstrated 
that the polarization correlations cannot be explained 
by any local realistic theory that uses hidden variables. 
Hence, together with other Bell test experiments, two- 
photon studies proved that nature indeed docs exhibit 
quantum-mechanical non-locality. These results con- 
tributed to the long-lasting historical debate of Einstein 
with Bohr and Schrodinger [l9|, [2(| who introduced the 
notion of entanglement for denoting nonproduct (pure) 
states. 



In the past, both experimental [17[ and theoretical 
[23l [HJ studies of 7 — 7 polarization correlation were 
mainly restricted to the 2s —> Is decay of neutral hy- 
drogen (or deuterium). Much less attention was paid to 
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two-photon transitions in other atomic or ionic species. 
With the recent advances in heavy-ion accelerator and 
trap facilities as well as in x-ray detection techniques, 
new possibilities arise to study spin-correlation phenom- 
ena in the decay of heavy, few-electron ions. In the 
medium- and high-Z domain, however, a proper anal- 
ysis of polarization quantum correlations requires de- 
tailed knowledge of relativistic and of many-electron ef- 
fects, as well as of those contributions that arise from the 
higher-order (non-dipole) terms in the expansion of the 
electron-photon interaction. 

In this contribution, we investigate quantum corre- 
lations between the polarization states of two photons 
emitted in the decay of few-electron ions. Most natu- 
rally, spin-correlation phenomena are described within 
the framework of density matrix theory. However, before 
we present details from this theory, we first summarize 
the geometry under which the two-photon decay is con- 
sidered in Section [H] In Section fill A[ then, the general 
expression for the spin-density matrix of the photon pair 
is derived, in terms of the initial populations of the ionic 
substatcs, and of the (second-order) transition ampli- 
tudes. The evaluation of these amplitudes in relativis- 
tic, second-order perturbation theory is thereafter dis- 
cussed for hydrogen- and helium-like ions. For the latter 
species, we make use of the independent particle model 
(IPM) which is appropriate for the analysis of bound- 
state transitions in the high-Z domain [U [22[. Apart 
from rigorous relativistic results, we also present simpli- 
fied expressions describing the photons' polarization state 
within the dipole approximation. These intuitive expres- 
sions, derived in Section IIIIB1 will enable us to under- 
stand the general behaviour of polarization correlations. 
In order to provide a quantitative description for these 
correlations, we briefly recall in Section IIVI the defini- 
tion of concurrence as measure of entanglement. Fully 
relativistic calculations of the concurrence are then per- 
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formed for the 2s\/2 — > lsi/2 transition in hydrogen-like, 
as well as ls 1/2 2s 1/2 -> ls^ 2 x 5 , ls 1/2 2s 1/2 3 S' 1 

^ s i/2 1 '^o an d lsi/2 2pi/2 3 -Po ~^ l s i/2 1 '5'o transitions in 
helium-like ions. Results of these calculations are dis- 
played in Section [VI and are compared to the predictions 
based on the dipole approximation. From this compar- 
ison we infer the twofold impact of relativity on polar- 
ization entanglement: apart from (i) the loss of purity of 
the photon states, (ii) the relativistic contraction of the 
wavcfunctions and the non-dipole contributions to the 
electron-photon interaction generally lead to the reduc- 
tion of concurrence; an effect that becomes prominent for 
heavy ions and high photon energies. A brief summary, 
together with some perspectives, is given in Section [VII 

Atomics units are used throughout the paper, unless 
differently stated. 



II. GEOMETRY OF THE SETUP AND THE 
PHOTON LABELING PROBLEM 

In order to analyze 7 — 7 polarization correlations, we 
first introduce the geometry of the two-photon emission. 
Since, for the decay of unpolarized ions, there is no direc- 
tion initially preferred for the overall system, we adopt 
the momentum of the "first" photon to coincide with 
the z axis which is also taken to be the quantization 
axis. Together with the direction of the "second" pho- 
ton, this axis defines the reaction plane (x-z plane). A 
single opening angle is required, therefore, to charac- 
terize the emission of the photons with respect to each 
other (cf. Fig. [J). 

Since the two photons are in a symmetrized state, it is 
a priori not possible to address them individually. How- 
ever, we can safely assume that photons observed by the 
detectors have definite energies and momenta, i.e. they 
collapse onto energy and momentum eigenstates. A clear 
identity can be given, therefore, to the photons [25| : the 
first (second) photon is that one detected by the detector 
A(B) (marked gray in Fig. [!} at a certain energy Wi( 2 ) 

and with momentum fci(2). By distinguishing in such a 
way the photons, we can use their polarization states in 
order to investigate the associated entanglement proper- 
ties. Indeed, such an analysis is possible since — in con- 
trast to the energy and momentum spaces — the photon 
spin state can be directly measured in any basis. 



III. THEORY 

A. Density matrix approach 

Having defined the geometry of the two-photon decay, 
we shall next recall the theoretical background needed 
to investigate the polarization of the emitted radiation. 
Most naturally, polarization-correlation studies can be 




FIG. 1: (Color online) General decay geometry. The z-axis 
is oriented along the direction of the photon with momentum 
fei, measured by the detector A. Together with the emission 
direction of the second photon, it defines the x — z reaction 
plane. The opening angle between both photons is denoted 
by 8, while the (linear) polarization angles of the photons, 
measured with respect to the ar-z-plane, are denoted by \i 
and X2> respectively. 



performed in terms of the system's density matrix. Since 
this approach was recently applied to describe two- 
photon transitions in hydrogen-like ions [23|, |24 |. here 
we may restrict ourselves to a short compilation of the 
basic formulas relevant for our further analysis. 

The initial state of the overall system, in our two- 
photon decay problem, is given by the photon vacuum 
I vac) = |0, 0)^ , and by the excited ion (or atom) in states 
I a,, Ji, Mi) with well-defined total angular momentum Jj 
and associated projection Mj onto the z-axis. Moreover, 
o,i is a collective label for all additional quantum num- 
bers required for a unique specification of the state. In 
particular, it characterizes electronic configurations that 
give rise to the state and, hence, provides its parity. Pi. 

The magnetic sublevel population of the ion in initial 
states is described as a statistical mixture, by the density 
operator 

Pi.ion = ^2 G M% \oti, Ji,Mi) (an, Ji, M;| , (1) 

Mi 

where Cm; denotes the population of the magnetic sub- 
state \o.i, Ji, Mi). Since in most (two-photon) experi- 
ments, the initially "prepared" excited ionic states are 
unpolarized, we fix the parameters as Cmi = l/(2Ji + 1). 
Such a realistic choice for the initial-state population has 
important consequences for the (spin) entanglement of 
emitted photon pairs. As we will see later, by introducing 
the incoherent mixture of initial magnetic substatcs ([1}, 
also the two-photon state's coherences are jeopardized 
and, hence, a loss of quantum correlations is induced. 

The density (statistical) operators of the initial and 
the final states of the overall system arc connected by 



3 



the standard relation I26H2 



pf = U pi,i on <g> p in , 



(2) 



where U is the evolution operator which accounts for 
the interaction of the ion with the radiation field. The 
final-state operator ^ describes both the de-excited 
ion in some state \ctf, Jf,Mf) and the two emitted pho- 
tons with momenta fci )2 and helicities Ai j2 . Owing to 
the transverse character of the electromagnetic radiation, 
these helicities, or projections of the photon momenta on 
their own directions of propagation, can take only two 
values Ai,2 = ±1. 

Instead of using the final-state density operator pf, 
it is often more convenient to work with its matrix rep- 
resentation, briefly referred to as the final-state density 
matrix. In the representation of the individual angular 
momenta this matrix reads: 



/; fciAi,fc 2 A 2 |/5/| /'; ki\' 1 ,k 2 X' 2 

ctfJfMf, fciAi,fc 2 A 2 \pf\atfJfMf', fciA' 1; fc 2 A 2 
1 



2Ji 



- dM k f t (Ai, M)M)t *(Ai, A 2 ) , (3) 

Mi 



where we employed Eq. ([1} in order to evaluate elements 
of the ionic initicil - stcitc matrix (o^t/jil^ jp^ion 
and introduced formal notation for the transition ampli- 
tude: 



M^ fe2 (Ax,A 2 ) 



J f M f ; fciAi,fc 2 A 2 U a, JiM i; 0,0) . (4) 



As seen from this expression, the M k f \ k2 (Ai, A 2 ) describes 
a transition between two bound ionic states accompanied 
by the simultaneous emission of two photons. 

The final-state matrix ([3]) still contains the complete 
information about the system and can be used to derive 
the properties of the photons or the residual ion. Assum- 
ing that the final magnetic sub-state of the ion remains 
unobserved in an experiment, we can derive the reduced 
density matrix p 7 which only describes the polarization 
state of the two photons, measured at a certain opening 
angle 6, with certain energies ui\ and u) 2 : 

(ki, Ai, fc 2 , A 2 | pf tl \k\, A' 1; fc 2 , A 2 ) 
= E (/; fc i' A i^2,A 2 |/5/ |/;fci,A / l! fc 2 ,A / 2 ) 



M, 



Af 



2.L 



Mi, Ms 



(5) 

where we introduced the factor N to ensure the proper 
normalization of the matrix, Tr(p 7 ) = 1. In what follows 
we will use this (reduced) matrix in order to analyze the 



polarization entanglement of the photons' pair. Before 
starting such analysis, we shall briefly discuss the compu- 
tation of the second-order amplitude (UJ) . Most naturally 
such amplitude can be evaluated within the framework 
of the second-order perturbation theory [|| [2(| [27| : 



A^ fc2 (Ai,A 2 ) = 

E v — Ei + ui 2 

(f\n(k 2 ,\ 2 )\is) H^(fci,Ai)|z) 

Ey — Ei + LU± 



(6) 



Here, \i) = \ai,J it Mi), \v) = \a v ,J v ,M v ) and |/) = 
\atf, J/,Mf) denote the solutions of Dirac's equation for 
the initial, intermediate and final ionic states, respec- 
tively, while Ei, E v and Ef are the corresponding en- 
ergies. Because of energy conservation, Ei and Ef are 
related to the energies lo\. 2 of the emitted photons by: 



Ei - Ef = wi + uj 2 



(7) 



From this relation, it is convenient to define the energy 
sharing parameter r\ = uji/(Ei — Ef), i.e. the fraction of 
the energy which is carried away by the "first" photon. 

In Eq. © TZ(k, A) is the transition operator that de- 
scribes the relativistic interaction of the electrons with 
the electromagnetic radiation. In velocity (Coulomb) 
gauge, this operator can be written as a sum of one- 
particle operators 



a m u\e 



(8) 



where a m denotes the vector of the Dirac matrices for the 
m-th particle, A\, m is the vector potential of the photon 
field and u\ is the unit polarization vector. For practical 
computations, it is convenient to decompose the vector 
potential A\, m into spherical tensors (i.e., into its electric 
and magnetic multipolc components). For the emission 
of the photon in the direction k = (9, 4>) with respect to 
the quantization z axis, such a decomposition reads: 



Ay 



^E E E* 

-L n p=0,l 



L '< [L 7 ] 1/2 {i\) p 



= 1 M, 



(9) 



where [L 7 ] = 2L 7 



1, k = \k\, D A j , is the Wigner 



rotation matrix of rank L 7 and the a p L ^ (k) refer to 

magnetic (p = 0) and electric (p = 1) multipoles, respec- 
tively. 

The great advantage of the multipole expansion ©, 
when comparing to the plane-wave formulation in the 
right-hand-side of Eq. (JSJ, is that it provides a radial- 
angular representation of the photon wavefunction. To- 
gether with the similar representations of the atomic 



4 



wavcfunctions it allows for significant simplification of 
the transition amplitude A^ fe2 (Ai,A 2 ) (see Ref. @ for 
further details). Moreover, Eq. §§§ gives a very useful 
tool for studying the multipole effects in the electron- 
photon interaction. If, for example, the summation in 
(|9]) is restricted to the term with L 1 = 1, p = 1, one 
obtaines the electric-dipole (El) contribution, while the 
Lj = 1, p = component provides the magnetic-dipolc 
(Ml) contribution, and so on. 



As seen from Eq. ([5]), the evaluation of second-order 
transition amplitudes requires the summation over the 
complete spectrum of the ion. Within the relativistic 
framework, such a computation is not a simple task 
since it includes a summation over the discrete part 
of the Dirac spectrum as well as an integration over 
the positive- and negative-energy continua. A number 
of methods have been developed over the past decades 
to compute Eq. ([5]) consistently. Apart from a direct 
summation over just few intermediate states which are 
close in energy to the states involved in the decay, the 
discrete-basis-set approach is widely employed nowadays 
in (relativistic as well as non-relativistic) second-order 
calculations. Within this approach, a finite set of dis- 
crete pseudostates is constructed from some basis func- 
tions and utilized for computing the transition amplitude 
-M.fi H, H(| ■ In the present work we use an alterna- 
tive, Green 's-function approach that helps to avoid the 
direct summation over the (virtual) intermediate states 
\v) = \a v ,J u ,M v ). In the framework of this alternative 
approach, moreover, we employ the Sturmian representa- 
tion [U of the radial components of the Green's function 
that allows the analytical evaluation of the transition am- 
plitudes ([5]) and, further down, of the entanglement mea- 
sures. 



In contrast to hydrogen-like ions, the relativistic 
second-order calculations for few-electron systems are 
more intricate, since one has to take into account 
electron-electron interaction effects. In the high-Z do- 
main, however, the radiative transitions in few-electron 
ions can be reasonably well understood within the IPM. 
This model, which takes the Pauli principle into account, 
is especially justified for heavy species, since the inter- 
electronic effects scale with 1/Z and, hence, are much 
weaker than the electron-nucleus interaction. The great 
advantage of IPM is that it allows the decomposition of 
the many-body, second-order transition amplitudes in 
terms of their one-electron analogs (see Ref. [22| for fur- 
ther details). In Section lVBl we will apply this approach 
for the computation of the spin-entanglement between 
photons emitted in the decay of heavy helium-like ions. 



B. Dipole approximation 

1. S — > S transitions 

The reduced density matrix ([5]) contains complete in- 
formation about the spin states of the photon pairs emit- 
ted in the decay of atoms or ions. Together with the tran- 
sition amplitude (|6]), it is suitable to explore two-photon 
transitions also in the high-Z domain, where relativistic 
and non-dipolc effects arc significant. However, before 
we perform such a fully-relativistic analysis, let us re- 
strict ourselves first to the non-relativistic dipole theory 
and derive approximate expressions for the description of 
the two-photon polarization states. As we will see later, 
this will provide intuitive insight into the entanglement 
properties of the photon pairs. Moreover, by comparing 
predictions of such a simplified dipole approach with the 
results of the fully-relativistic theory, we will be able to 
identify the relativistic and multipole effects in the two- 
photon transitions. 

By making use of the non-relativistic dipole approx- 
imation for the electron-photon interaction and by re- 
stricting the intermediate-state summation to states 
\v) = \a v , J Vl M v ) with definite momentum J„ and par- 
ity P„, defined by dipole selection rules, it is possible to 
express the two-photon density matrix ([5]) in the form 



(fci,Ai,fc 2 A 2 |/5 7 \ki, Ai,fc 2 , A' 2 ) 
«CA 1 A 2 A' 1 A 2 D^ 2 (x'y'z' ^xyz)x 

x (1, Ai, 1, -X[\L, /n) (1, -A 2 , 1, X' 2 \L, fi 2 ) x 





Jf 






Jv 






1 





1 1 L 



1 1 L 



(10) 

where we employ the standard notation for the Wigner 
6j- and 9j-symbols [HI, and C is a normalization con- 
stant that absorbs the radial parts of the (dipole) transi- 
tion amplitudes Ai f{. The final-state density matrix (fT0|) 
depends, therefore, only on the symmetry of the initial 
and final ionic states as well as on the photons' helicities. 



As mentioned above, Eq. (fTOf can be applied for the 
analysis of only those transitions that proceed — within 
the non-relativistic picture — via intermediate (virtual) 
states \a„, J Vl M v ) having one particular value of to- 
tal angular momentum J„ and parity P„. This is the 
case of the lsi/ 2 2sj/ 2 1 Sq — > ls^/ 2 transition in 
helium-like ions, for which the intermediate-state sum- 
mation in the amplitude (j6|) is restricted to the n v 1 P\ 
levels only, if one treats electron-photon interaction in 
the dipole approximation. Our approach is also justi- 
fied for the 2si/ 2 — > lsi/2 decay in hydrogen-like ions, 
since the intermediate states n v p\/2 and n„p 3 / 2 are de- 
generate in the non-relativistic limit. In the following 
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FIG. 2: (Color online) Concurrence of two photons emitted in the 2si/2 — > lsi/2 decay of neutral hydrogen and hydrogen-like 
xenon and uranium ions. Results of the non-relativistic dipole approximation (dashed line) and the rigorous relativistic theory 
(solid line) are shown, for two relative photon energies r) = 1/16 (upper panel) and <t] = 1/2 (lower panel). 



analysis, therefore, we shall restrict the discussion of the 
non-relativistic dipole approximation (|10[) to these two 
transitions. 



By inspecting Eq. (fTU|) for the cases of the (non- 
relativistic) 2s 1 / 2 —> lsi/2 an d ^Si/2^ s i/2 So 
l_ s i/2 1 Sq transitions, it can be proven that the density 
matrix p 7 fulfills the relation Tr(pi) = (Tr(/S 7 )) 2 = 1 
and, hence, represents a pure quantum-mechanical spin 
state of the emitted photons. This pure state can be 
described by the state vector: 



1 



2Vl + cos 2 9 
+(cos0 + l)(|+-> + |-+» 



(cos0 -l)(|++) + | )) 



(11) 



as directly derived from Eq. (fTU|) . In this expression, 
the pre-factor arises due to the normalization condition, 
(vp | = 1 and, for the sake of brevity, we use the no- 
tation |±) = |A = ±1). Here and henceforth, whenever 
a state vector describing both photons appears, the first 
(second) index has to be attributed to the first (second) 
photon; the photons being identified and detected ac- 
cording to Sec. HIl For the particular case of back-to- 
back photon emission (9 = tt), the vector (fTTj) further 
simplifies to a Bell state: 



l*> = i(l++> + |— » 



(12) 



As seen from this expression, for the opening angle 9 = tt 
photons can only be detected having the same helicity. 
In the past, such a quantum correlation between the 
photons emitted in the 2s]y 2 — > lsi/2 decay of atomic 



deuterium has been employed for verifying a violation of 
Bell's inequality QJ]. 

Instead of the helicity basis |A = ±1), it is often more 
convenient to analyze the polarization correlations of two 
photons in terms of their linear polarization unit vectors. 
These vectors are defined in the plane perpendicular to 
the photon propagation axis and can be obtained by the 
standard transformations 1321: 



l*> = 75 (l+> + |-» > 

\y) = is(- !+> + !-», 



(13) 



and can be used to re-write the state vector (jTTJ) in the 
xy-representation as: 



I*) 



1 



Vl + cos 2 9 



\yy) + cos 9 \xx) 



(14) 



Since we just change the basis, the opening angle 9 = tt, 
again corresponds to a Bell, i.e. maximally entangled, 
state. In contrast, in this representation (|14[) we imme- 
diately see that perpendicular emission under 9 = tt/2 
results in a product (non-entangled, or separable) pho- 
ton spin state. 

The quantitative analysis of entanglement for the emit- 
ted photon pair will be performed in the following Sec- 
tions, based on our general expression (|10p . as well as 
on the dipole approximation ([l"2"j) - (fl"4"|) . Before we start 
with such analysis, let us still discuss some basic prop- 
erties of the spin-state (JUJ). In this way we recall the 
typical "experimental scenario" in which both photons 
arc detected by polarimctcrs whose transmission axes are 
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characterized by the angles xi an d X2 with respect to the 
reaction (x-z) plane (see Fig. [TJ. Eq. (fT4"|) predicts that 
after the first photon has been detected by the detector 
A (the first detector) with some defined (linear) polar- 
ization angle xij the second photon collapses onto the 
vector: 

I*) ->• 7V[sinxi \y) + cos6»cosxi \x) ] , (15) 

with N some normalization factor. It follows from 
Eq. (Ti~5| that the second photon is then found in a lin- 
early polarized state. The direction of this linear polar- 
ization, characterized by the angle X2, depends on the 
opening angle 6 and on the polarization angle xi of the 
first photon: 

tanx2 = —^-7 tan xx • (16) 
cos# 

As we will show later, such a definite (except for the 
opening angle 9 = tt/2) correspondence between the lin- 
ear polarizations does not generally imply maximal en- 
tanglement of the photon pairs. To understand this issue 
and to quantify the degree of entanglement we shall in- 
troduce, in Section [IV! the concurrence measure. 

2. P — > S transitions 

In contrast to the S — s- S transitions from above, the 
non-relativistic dipole approximation (|10[) - (|14[) cannot 
be applied for the analysis of the lsx/2 2px/2 3 -fo — ► 
l' s i/2 ^° decay of helium-like ions. The principal rea- 
son for this failure is that the leading (electric-magnetic 
dipole) E1M1-M1E1 ls 1/2 2p 1/2 3 P -> ls\ /2 1 S transi- 
tion may proceed either via intermediate lsx/2 n u s i/2 3 Si 
or ls 1 / 2 n u p 3 Pi states, thus giving rise to a "double- 
slit" picture. By taking into account such a Young-type 
interference and by restricting ourselves to the dipole 
(E1M1-M1E1) terms in the electron-photon interaction, 
we again find the photon pair in a pure state: 

I*} = c(-E( 7 y)sin 2 ^|++)+A(7 ? )cos 2 ^|-+) 
-A(7?) cos 2 f|+-> + £(??) S in 2 § I— >)■ 

(17) 

Here, C is the normalization constant, and the energy- 
dependent functions £(?/) = S E imi{^i) + S E imi{u 2 ) + 
Smiei(wi) + SAfiBi(u> 2 ) and A (77) = Seimi(ui) - 
Seimi(u2) — Smiei(ui) +Smiei(u2) are given in terms 
of the multipole, second-order reduced transition ampli- 
tudes Sl iPi ,l 2 p 2 ( U} ) ( see R- ei - [Hi f° r further details). 

As seen from Eq. (|17|). the spin-state of the photons 
emitted in the lsx/2 2px/2 3 Po — > ls 2 / 2 So transition 
depends on the energy sharing 77. No simple analytical 
expression for this dependence can be derived in general, 
owing to the complicated structure of the functions A (77) 
and £(??). However, if both photons carry away the same 



fraction of the energy, u>i — us 2 , the function A (77 = 0.5) 
vanishes, and the vector (j!7[) represents a maximally en- 
tangled (Bell) state: 

I*) = ^= ( - + I— > ) = ( \xy) + \yx) ) . 

(18) 

By comparing this expression with Eq. (|14[) . one can see 
that polarization properties of Pq — > Sq and So — > So 
(as well as non-relativistic Si/ 2 — > Sx/2) transitions are 
rather different: while the photons emitted in the S — > S 
transitions can be detected having parallel linear polar- 
ization vectors, the P — > S decay should result in emis- 
sion of the photon pair with orthogonal linear polariza- 
tions. Moreover, no angular dependence arises in the 
state vector (|18l) implying maximal entanglement be- 
tween the photons' spins, irrespective of the particular 
decay geometry. 



IV. ENTANGLEMENT OF THE TWO-PHOTON 
STATE 

We are ready now to discuss the concept of entangle- 
ment for the emitted photon pairs and to introduce a 
proper measure for it. In this way, let us first return 
to the full two-photon state which accounts not only 
for the spin but also for the spatial degrees of freedom. 
Within the non-relativistic dipole approximation (|14l) . 
such a state before its detection reads: 

m = nJ dujiduj 2 5(uji + uj 2 - A u )/(wx) 

x f d0xd0 2 (\yy) 

+ 008(61! - e 2 ) \xx) )+(1h2), (19) 

where (1 «-> 2) denotes the previous terms but with all 
particles' labels exchanged and the state is normalized by 
virtue of the constant N. Moreover, f(u>i i2 ) is the energy 
probability density function of the decay, A w = Ei — Ej 
is the transition energy and 6^2 are the angles which 
address the position of the first and second photon in 
the reaction plane, respectively. Due to the integral over 
angles and energies, the above state can be written as 
product state neither in the energies nor in the emission 
angles or in the polarization. It can hence be seen as 
highly entangled, in general. 

In order to rigorously discuss entanglement, due to the 
identity of particles, a degree of freedom for discrimi- 
nation is needed. The energy of the photons and their 
opening angles are an appropriate choice, since one nat- 
urally projects onto energy and momentum eigenstates 
in the experiment. In the coincidence experiment dis- 
played in Fig. [I] the two-photon state collapses onto a 
state with definite momenta, while the polarization can 
be measured in any basis. If the energies of the photons 
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are equal and the emission directions exactly the same, 
we are not able to identify two separated particles be- 
tween which entanglement may be defined. As long as 
this is not the case, we identify the particle projected on 
the two energies and angles as the two distinct entities 
to which we can safely assign an entanglement measure 
[25j . Hence, even though we start with a rather complex 
state of identical particles in which no physical subsys- 
tem structure is apparent, we can effectively deal with 
the two-qubit system of polarized photons projected on 
energy and momentum states. 

Having clarified the concept of the two-photon entan- 
glement, we shall introduce now its quantitative measure. 
For a photon pair, that can be seen as a "two-qubit" 
system, it is very convenient to describe the degree of 
entanglement by means of the Wootter's concurrence C 
[33[ | . For an arbitrary two-qubit system described by the 
density operator p the concurrence is defined as 



U 



C 



0, 



(20) 



where J~el are the square roots of the eigenvalues of the 



operator p(&rp <g> a^')^^^' ® m descending or- 
der, p* is the complex conjugate of p, and cr^ 1 ' 2 ^ = S'y" 2 ^ 
are the Pauli matrices acting on the first and the sec- 
ond qubit, respectively. Before we discuss further the 
properties of the concurrence C, let us first recall that it 
quantifies correlations that can be fully attributed to the 
entanglement. Bi-particle states with vanishing concur- 
rence can still exhibit correlations which are, however, 
not of quantum nature. 

Definition (f2"0"|) can be simplified further if applied to a 
pure quantum-mechanical state described by a ket vector 



r (2) 



= C aa \aa) + Cab \ab) + 
C ba \ba) + C bb I bb) 



(21) 



where a, b are arbitrary two-dimensional basis states and 
Cij are complex numbers. For this state, the concurrence 
reads 



C — 2 C aa Cbb ~ C a bC\ 



bu 



(22) 



By using this expression and Eq. 
obtain the analytical expression 



C(6) 



cos ( 



1 + cos 2 



T4|) . we immediately 



(23) 



for the spin-entanglement of the photons emitted in the 
2s 1 / 2 I.S1/2 and ls 1 / 2 2s 1/2 1 S -> ls| / 2 ^0 tran- 
sitions. We remind that Eq. ([23]) is obtained within 
the non-relativistic dipole approximation and should be 
questioned in the high-Z domain, where higher-order 
and relativistic effects can play a significant role. To 
explore the influence of these effects on the photon spin- 
entanglement, we will compare in the next section the 
predictions obtained from Eq. (|23[) with the rigorous rel- 
ativistic calculations based on Eqs. (O-© and (|20p. 
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FIG. 3: (Color online) Purity 1)25|) of the two-photon state in 
the 2s 1 / 2 — ► lsi/2 decay of hydrogen-like uranium. Results of 
the non-relativistic dipole approximation (dashed line) and 
a rigorous relativistic treatment (solid line) are shown, for 
relative photon energies rj = 1/16 (upper panel) and 77 = 1/2 
(lower panel). 



V. RESULTS AND DISCUSSION 
A. Hydrogen— like ions 

1. Polarization entanglement 

After discussing the theoretical background of two- 
photon polarization studies, we are prepared now to ana- 
lyze the influence of the relativistic and higher-multipolc 
effects on quantum correlations between the emitted par- 
ticles. We start our analysis with the 2si/ 2 — > lsi/2 decay 
of hydrogen-like ions, that is well established both in the- 
ory and in experiment. As shown above, the polarization 
properties of this transition can be described — within the 
non-relativistic dipole approximation — by the state vec- 
tor ([TJI and, hence, by the degree of entanglement (|23|) . 
The theoretical predictions, obtained within such a non- 
relativistic approach, are displayed in Fig. [2] for the de- 
cay of neutral hydrogen as well as hydrogen-like xenon 
Xe 53+ and uranium U 91+ ions and are compared with the 
results of the rigorous relativistic treatment. For the lat- 
ter, one deals with the relativistic Dirac's wave-functions 
and includes, in addition, the full interaction between the 
electron and the radiation field in the amplitude (0 . Rel- 



8 



Xe 



52+ 



u 



90+ 




60 



120 



180 




60 



120 




0.8 



0.6 



0.4 



0.2 



180 



1 1 1 1 1 










A / 






a // 






\\ // 






V/ 





60 120 180 60 120 180 

Opening angle 9 (deg) Opening angle 9 (deg) 



FIG. 4: (Color online) Concurrence of two photons emitted in 
the ls 1 /22p 1 / 2 S Po — > lsi/2 1 5'o decay of helium-like xenon 
and uranium ions. Results of the E1M1 dipole approximation 
(dash-dotted line) and of rigorous relativistic theory (solid 
line) are shown, for two relative photon energies r\ — 1/10 
(upper panel) and r\ = 1/4 (lower panel). 



maximally entangled Bell states |*) = (|H — ) + |H — ))/V2 

for 9 = 0, and |*) = (|++) + | ))/ v / 2 for 9 = tt, as 

predicted by Eqs. (IT]) and (|T2|) . 

Similar to the cases of parallel and back-to-back pho- 
ton emission, the conservation condition Ai = — M 72 with 
the helicity of the first photon being Ai = ±1 may help to 
understand the behaviour of the entanglement measure 
C(9) at 9 = ir/2. This will require us to return to Eq. © 
which simplifies — within the dipole approximation — to: 



X x = -Vfo J2 Mm>)4/ 7 a(^ (24) 

M T = -1 

where d} M x (9) is the Wigner's (small) <i-matrix whose 
properties are discussed in detail in [28J. For the open- 
ing angle 9 = tt/2, the elements of this matrix are: 
d\i = d\_ 1 = dl n = d 1 _ 1 _ 1 = 1/2, implying, together 
with Eq. (|6]) and with the fact that the ionic states are 
spherically symmetric, that the probability for the sec- 
ond photon to have projection M 72 = =pl on the quan- 
tization axis of the overall system is independent of its 
helicity A2. Thus, no correlations between the polariza- 
tion (spin) states of the emitted photons appear for the 
perpendicular emission, leading to the vanishing entan- 
glement C(tt/2) = as displayed in Fig. [2] 



Purity of the two-photon state and impact on 
entanglement 



ativistic as well as non-relativistic dipole calculations of 
the concurrence (|2"0|) - (f2"2"]) are performed at two relative 
photon energies i] = 1/16 (upper panel) and 77=1/2 (lower 
panel). As seen from the figure, in case of equal energy 
sharing (7/ = 1/2), both approaches yield almost identical 
results along the entire isoelectronic sequence. Our cal- 
culations show that, while being maximal for the parallel 
(8 = 0) and back-to-back (9 = tt) photon emission, the 
concurrence vanishes at the opening angle 9 = tt/2. This 
behaviour is well understood from Eqs. (|lip and (fTl)) 
as well as from the conservation of the projection M to t 
of the total angular momentum J t ot of the overall sys- 
tem "ion + two photons" . Namely, if no electron-spin 
flip were to occur during the 2s l / 2 — > ls 1 / 2 decay and 
assuming nuclear spin, the conservation law enforces 
that the change of the projection of the ion's total angu- 
lar momentum relative to the quantization axis z (chosen 
along the momentum of the first photon) would be given 
by Mi — Alf = = Ai + Af 72 . In this expression, Ai is 
the helicity of the first photon and Af 72 is the projection 
of the angular momentum of the second one. For pho- 
tons emitted in parallel or back-to-back, this projection 
is M 72 = A2 and M 72 = — A2, correspondingly, thus lead- 
ing to the conditions Ai = — A2 or Ai = A2. Moreover, 
owing to the spherical symmetry of s-ionic states there 
is an equal probability of emission of the "first" photon 
with helicity Ai = +1 or —1. This immediately implies 



As seen from the top panel of Fig. O the accuracy of 
the non-relativistic approximation (|23p becomes gener- 
ally worse if one of the photons has a significantly larger 
energy than the other one. For the 2s i / 2 —> lsi/2 tran- 
sition of hydrogen-like uranium, for example, the non- 
relativistic dipole approximation overestimates the con- 
currence measure by about 10% for forward as well as 
backward opening angles, for an energy sharing ?/ = 1/16. 
In order to understand better such an energy-dependent 
behavior, we study the purity of the two-photon polar- 
ization state, defined as 

7> = (25) 

where p 1 represents the photon density matrix (|10p . 
The purity varies from (completely mixed state) to 1 
(pure state). In Fig. [3] we display the purity V for the 
2si/2 l s i/2 decay of U 91+ for two relative photon en- 
ergies: 77 = 1/16 (upper panel) and 1/2 (lower panel). As 
seen from the figure, the purity strongly depends on the 
energy sharing parameter: while the purity of the two- 
photon state is always > 0.987 for an equal energy sharing 
?; = 0.5, it is significantly reduced for 77 = 1/16. The loss 
of purity can be attributed to the spin-orbit coupling in 
hydrogen-like ions as well as to the magnetic terms in 
electron-photon interaction. Both these relativistic ef- 
fects increase with the nuclear charge Z and with the 
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FIG. 5: (Color online) Concurrence of two photons emitted 
in the lsi/2 2si/2 3 Si — > lsf/ 2 1 So decay of helium-like xenon 
and uranium ions. Results of the E1E1 dipole approximation 
(dash-dotted line) and the rigorous relativistic theory (solid 
line) are shown, for two relative photon energies r\ = 1/10 
(upper panel) and r\ = 1/4 (lower panel). 



photon energy u>. They lead to the fact that the decay 
of the unpolarized and, hence, mixed 2si/ 2 level results 
in the emission of photons characterized by a partially 
mixed state. Due to complementarity of entanglement 
and mixedness/impurity [291 ] . such a loss of purity causes 
the reduction of the concurrence measure that can be 
observed in the top panel of Fig. [2j Despite such a re- 
duction, there are still quantum correlations between the 
polarization states of the photons. 



Helium— like ions 



1. lSi/2 2S!/ 2 So 



So transition 



In contrast to the 2si/2 — > lsi/2 decay of one-electron 
systems, the ls 1 / 2 2s 1 /2 1 S'o — > l s i/2 1 So transition 
in helium-like ions always proceeds between the pure, 
J = 0, quantum-mechanical states (cf. Eq. ([1])). There- 
fore, the spin-state of the photon pair emitted in such 
a transition cannot be but pure at any energy shar- 
ing. The concurrence measure C{9) of such a pure 
state calculated within the exact relativistic theory turns 
out to be almost identical to the dipole approximation 
(|2U|) . The deviation between both predictions does not 
exceed 1%, even for the heaviest helium-like ions and 



arises due to the higher, non-dipole terms in electron- 
photon interaction ((HJ. By comparing this prediction 
with the calculations performed in the previous section 
for hydrogen-like ions, we again argue that the reduction 
of the spin-entanglement between the photons emitted in 
the 2s 1 /2 — > lsi/2 transition shall be mainly attributed 



to the loss of purity of ionic states. 



2. lsi/2 2pi/2 3 -Po — > ls?/2 "0 transition 



After the above, short discussion of the 
lsi/2 2si/2 So — > l s r/2 transition, we 

turn now to explore quantum correlations in the 
ls 1 / 2 2p 1 / 2 3 -Fb l s i/2 1, S'o decay. As one can expect 
from the spin-state vector (JTTJ) , derived in leading order, 
"electric and magnetic" dipole approximation, these 
correlations should differ from those predicted for the 
S — > S cases. Indeed, by inserting the vector (fT?|) 
into Eq. l[2"g]). we obtain — within the dipole (E1M1) 
approximation — the concurrence measure as: 



C(e,r,) = 



A 2 (77) cos 4 f -S 2 (r/) sin 4 f 



A 2 (77) cos 4 ! + S 2 (77) sin 4 ! 



(26) 



In contrast to the polarization entanglement (|23p be- 
tween the photons emitted in the S — > S transitions, 
the concurrence turns out here to depend on the pho- 
tons' energy sharing. To better understand such a de- 
pendence, we display the entanglement of the photons 
emitted in the lsi/2 2pi/ 2 3 fo —> l- s i/2 1 ^o decay of 
helium-like xenon and uranium ions in Fig. [?] Calcu- 
lations arc performed for two relative energies: r/ = 1/10 
and 1/4 both within the dipole approximation (|26[) and 
by using the exact theory which accounts for the higher 
multipole channels. As seen from the figure, both the- 
oretical approximations predict maximal entanglement, 
C = 1 , for the parallel and back-to-back photon emission 
at any energy sharing 77; a feature that could be expected 
from the conservation laws. On the contrary, the "criti- 
cal" opening angle 9 C at which the concurrence vanishes, 
C(0 C ) = 0, varies with the relative photon energy. By 
inspecting Eq. ([26| we find the following relation for this 
angle: 



tan 



|A(t?)| 

\nv)\ 



(27) 



It follows from this expression that for any non-zero val- 
ues of the functions A (77) and £(77) there exists one single 
critical angle O c (rj), as can also be seen, for example, from 
Fig.Hl 

If the photons, emitted in the ls 1 / 2 2p 1 / 2 3 Po — ^ 
I s 1/2 decay, carry away the same portion of en- 
ergy, u)\ = uj 2 , the function A (77) turns out to be zero 
and Eq. (|27|) cannot be applied for the determination 
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of the critical angle 9 C . As can be seen from Eq. (fT8|) 
and Eq. (|26p . in this case the photons' state is max- 
imally entangled (Bell state) for any opening angle, 
i.e. C(9, 0.5) = 1. This behaviour differs from that of 
lsi/ 2 2s 1 / 2 So — > ls^ 2 1 5'o as well as 2si/ 2 — > lsi/ 2 
transitions for which no correlations appear at the open- 
ing angle 9 ~ ir/2. In order to understand the reason for 
this difference, we shall return to Eqs. ([B])-©- By mak- 
ing use of these expressions and of the properties of the 
Wigner matrices we re-write the two-photon transition 
amplitude in terms of reduced matrix elements as: 

7wff 2 (A l7 A 2 )cx (A 1 +A 2 )(W/iM + 5 M i£iM) , 

(28) 

where, for the case of equal energy sharing, ui\ = lo 2 = w. 
It follows from Eq. (|2"8"|) that apart from the conservation 
of the projection of the total angular momentum Jtot, dis- 
cussed in Section IV Al an additional selection rule arises 
for the ls 1 / 2 2p 1 / 2 3 Po — > l s i/ 2 that forbids emis- 
sion of the photons with opposite helicities. Together 
with the equal probabilities of the spin-states |++) and 

| ), this selection rule implies the Bell state (TT5|) and, 

hence, maximal engagement of the photons' state. 



3. Incoherent preparation of the ions: 
2si/2 3 Si lsi/2 1 transition 

Until now our discussion of the two-photon decay of 
helium-like ions was restricted toJ = 0— >J = transi- 
tions. In this case, both initial and final ionic states are 
pure along the entire isoelectronic sequence, and, con- 
sequently, the two-photon states are pure as well. In 
order to underline again the effect of the loss of pu- 
rity on the quantum correlations, we study the two- 
photon decay of the unpolarized ls 1 / 2 2s 1 / 2 Si state. 
In Fig. [5j we display the degree of spin-entanglement for 
the ls 1 / 2 2s 1 / 2 3 Si — > ls^ 2 So transition in helium- 
like xenon and uranium ions. Again, the exact rela- 
tivistic calculations are compared with the predictions 
of electric dipole (E1E1) approach for two relative en- 
ergies rj = 1/10 and 1/4. As seen from the figure, the 
general behaviour of the measure C3 5 _j.i5 is very sim- 
ilar to that of the ls 1 / 2 2s 1 / 2 1 So — > ls^ 1 Sq an d 
2si/ 2 — > lsi/ 2 transitions. Namely, the concurrence 
changes from C3g ^ig = 1 for the parallel photon emis- 
sion down to zero at 9 = tt/2 and back to a maximum 
entanglement for 9 = ir. Similar to the discussion in 
Section IV Al this can be easily understood if one ap- 
plies again the momentum projection selection rules and 
Eq. ([M]). In contrast to the lsi/ 2 2s 1 / 2 1 S -> ls^ 2 1 S 
and 2si/2 — > lsi/ 2 transitions, however, the degree of 
entanglement C3 Sa ^i Sa drops down much faster for the 
forward < 9 < tt/3 and backward 2ir/3 < 9 < ir angles; 
an effect that can be understood if we remember that the 
initial ionic state is prepared in an unpolarized (mixed) 
state. 



As one can see from Fig. [5j spin entanglement for the 
ls 1 / 2 2s 1 / 2 3 S\ — > ls^ 2 1 So transition is very sensitive 
to higher multipoles in the electron-photon interaction. 
This is a direct consequence of a strong suppression of 
the E1E1 decay channel caused by the symmetry prop- 
erties of the multi- pho ton systems as described by Bose 
statistics (see Refs. |34H37| for further details). The non- 
dipole contributions become more significant for heavier 
ions and with increasing energy sharing rj. < n < 0.5, 
and result in an asymmetric shift in the concurrence. 



VI. SUMMARY 

In summary, the two-photon decay of few-electron 
ions has been investigated within the framework of den- 
sity matrix and second-order perturbation theory. Spe- 
cial attention in our study has been paid to the quan- 
tum correlations between the spin states of the emitted 
photons. By making use of the non-relativistic dipole 
model, we derived a simple analytical expression for such 
spin-entanglement if observed in 2si/ 2 — > lsi/2 an d 
ls 1 / 2 2s 1 / 2 1 So —> ^- s \/ 2 Sq transitions in hydrogen and 
helium-like ions, respectively. By comparing predictions 
of the dipole approximation with the fully-relativistic 
calculations, we were able to explore the influence of the 
relativistic effects on the photon polarization properties. 
In particular, we observed a reduction of entanglement 
that becomes more sizable for high-Z systems and that 
can be attributed to the loss of purity of the two-photon 
spin states induced by higher multipole contributions and 
the relativistic contraction of the wave-functions. 

Beside the well-established 2s!/ 2 — s- lsi/ 2 and 
lsi/ 2 2s!/ 2 1 So — > ls^ 2 1 So transitions, entanglement 
studies were also performed for the ls 1/ / 2 2s 1 / 2 3 Si — > 
i s i/2 1 S'o and lsi/ 2 2p 1 / 2 3 P — s- ls^ 2 x Sq decays of in- 
termediate and high Z helium like ions. Based on the in- 
dependent particle model, which is a good approximation 
for the analysis of the bound-bound transitions in heavy 
atomic systems, we found that the concurrence is very 
sensitive to the relative energy of emitted photons, as well 
as to the higher multipole contributions to the electron- 
photon interaction. The strongest non-dipole effects 
have been identified for the ls 1 / 2 2s 1 / 2 3 Si —> lsy 2 1 So 
two-photon transition for which the E\E\ decay channel 
is forbidden due to symmetry properties of the system. 

Our theoretical analysis, performed for the intermedi- 
ate and high Z domain, underlines the importance of 
detailed knowledge on the electronic structure of ions 
(atoms) for a better understanding of two-photon entan- 
glement. This is rather different from the earlier studies 
on the decay of light neutral atoms fl6l - fT8T ] where the 
quantum correlations between the photons could be pre- 
dicted solely from angular-momentum conservation. In 
contrast, for heavy ions, entanglement properties of the 
photon pairs are governed by the complicated interplay 
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between these conservation laws and relativistic as well 
as many-body effects. 

The spin quantum correlation studies reported in the 
present work will help to understand the outcome of 
the future 7 — 7 coincidence measurements. Owing to 
the recent developments in x-ray polarization detectors, 
these measurements are likely to be performed at the 
GSI storage ring in Darmstadt. Apart from the analy- 
sis of relativistic and quantum electrodynamics (QED) 
effects, parity-violation (PV) corrections to the photon 
spin-entanglement can be observed in such coincidence 
experiments. A theoretical analysis of these PV phenom- 



ena is currently under way. 
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